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C. D.U. 624.134,25: 624.1431.36 
MÉTODO PARA AMOSTRAGEM INTEGRAL DE MACIÇOS 
ROCHOSOS 
por MANUEL ROCHA 
Director, Laboratório Nacional 
de Engenharia Civil, Lisboa 
RESUMO SYNOPSIS 


E apresentado novo método de amostragem com 
o qual é possivel obter tarolos que representam a totali- 
dade dos maciços rochosos, incluidas as fracturas de 
qualquer natureza e os seus enchimentos, assim como 
sonas de muito fraca resistência. Além disso, os tarolos 
obtidos são orientados. 

Us tarolos permitem a determinação das proprieda- 
des do maciço rochoso, em particular das fracturas 
e enchimentos. 


1. INTRODUÇÃO 


Apesar de já há muito tempo ser reconhecido 
o papel primordial desempenhado pelas diacla- 
ses, falhas e outras fracturas no comportamento 
dos maciços rochosos, o facto é que até hoje não 
foi possível, que saibamos, desenvolver técnica 
de amostragem em furos de sondagem que per- 
mitisse obter informação satisfatória sobre as 
características dessas discontinuidades, em par- 
ticular dos seus enchimentos. De um modo ge- 
ral, a amostragem que se consegue não é repre- 
sentativa das zonas de menor resistência, afinal 
aquelas que mais influenciam a deformabilidade 
e a resistência dos maciços rochosos. Na pros- 
pecção mineira o conhecimento das fracturas 
e seus enchimentos e de zonas de fraca resistên- 
cia é também muitas vezes de importância 
capital. 

É situação corrente serem despendidas somas 
avultadas em trabalho de sondagem sem se con- 
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À new sampling method is presented which yields 
core samples representative of the whole rock mass, 
including fractures of any nature and their infillings, 
and also of the zones with very low strength. Addi- 
tionally the core samples are oriented. 

By means of the core samples it is possible to deter- 
mine the properties of the rock mass, particularly of 
Jfractures and infillings. 


seguir obter as informações necessárias, o que 
cria muitas vezes verdadeiro sentimento de frus- 
tração. No caso de estar em jogo o comporta- 
mento mecânico, como os tarolos obtidos são 
representativos das zonas de melhor qualidade, 
as sondagens conduzem mesmo por vezes a um 
julgamento optimista das características do ma- 
ciço, o que pode acarretar as mais graves con- 
sequências. 

As dificuldades que se acabam de apontar 
levaram ao desenvolvimento de diversas técnicas 
de observação das paredes dos furos de sonda- 
gem, por meio de periscópios e câmaras fotográ- 
ficas e de televisão. Essas técnicas apresentam 
contudo severas limitações quer quanto à natu- 
reza e qualidade da informação obtida quer às 
condições de utilização. 

Em face da presente situação fomos conduzidos 
a procurar um método de amostragem que per- 
mitisse obter amostras intactas integrais, isto é, 
representativas da totalidade do maciço rochoso, 
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o que na nomenclatura habitual se pode exprimir 
por uma percentagem de recuperação de 100 */o. 
O método desenvolvido fornece ainda a orienta- 
ção das amostras. 


2—- NOVO MÉTODO DE AMOSTRAGEM 
DESENVOLVIDO 


O novo método de amostragem consiste essen- 
cialmente em obter tarolo do maciço rochoso 
previamente armado por meio de varão o qual 
assegura a solidariedade e a integridade do mate- 
rial extraído. O método poderá ser aplicado ao 
longo de todo um furo de sondagem ou sômente 
nas zonas onde as condições do maciço ou as 
exigências da obra a construir aconselhem a sua 
utilização. 

As operações a efectuar para a obtenção de 
amostras integrais orientadas são as seguintes: 


a) Abertura do furo de sondagem (Fig. 1 a). É aberto 
furo de sondagem, com a inclinação desejada, 
até à profundidade a partir da qual se pretende 


| PEÇA DE“ 
LIGAÇÃO 


he E + ar | 
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Extarosesrenk 


A COAXIAL 
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Fig. 1 


fazer a amostragem integral. O diâmetro D do 
furo será o adequado ao diâmetro do tarolo que 
se deseja obter. Temos sobretudo adoptado o diá- 
metro de 7,5 cm, comum nos trabalhos de pros- 
pecção, mas já se foi até 17 cm. 
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b) Abertwa do furo para alojamento do varão 
(Fig. 1b). No fundo do furo anteriormente refe- 
rido é aberto furo coaxial com diâmetro de cóm- 
primento | igual ao do tarolo a obter. Convém 
que o diâmetro d seja pequeno com vista a 
aumentar a espessura da amostra armada, mas 
por outro lado é necessário que o furo permita 
alojar varão com rigidez que evite deformações 
inconvenientes do tarolo no decorrer das opera- 
ções, adiante referidas, de sobrecarotagem, de 
extracção do amostrador e de colocação nas 
caixas de acondicionamento. Para furos de son- 
dagem com D==7,5cm tem sido usado com 
maior frequência o valor d ==2,6 cm, o que con- 
duz a tarolos com uma espessura de cerca de 
1,7 cm. Quanto ao comprimento | do furo, con- 
vém que seja grande a fim de reduzir o número 
de operações de amostragem para um mesmo 
comprimento total, e portanto o custo do traba- 
lho. O valor de | é por outro lado limitado pelo 
risco de, na operação de sobrecarotagem, se veri- 
ficar desvio exagerado na coaxialidade dos dois 
furos, podendo mesmo vir a ser interceptado 
o varão de reforço. Temos até agora obtido 
sobretudo tarolos com o comprimento de 1,5 m 
mas não prevemos que, na maior parte dos casos, 
haja dificuldade em ir bastante além deste valor. 

O furo para alojamento do varão pode ser 
aberto com coroa de diamante de tipo habitual 
ou com coroa de secção cheia. O primeiro tipo 
tem a vantagem de se obterem tarolos e portanto 
uma primeira informação sobre o maciço rochoso 
que pode ter utilidade para a condução das ope- 
rações subsequentes de obtenção da amostra 
integral. Por outro lado, usando coroas de secção 
cheia está-se menos sujeito ao risco de a coroa 
se encravar e de ao ser retirada ficar quantidade 
apreciável de detritos no fundo do furo. 

A coaxilidade do furo em relação ao furo de 
maior diâmetro é obtida por meio de guia que se 
apoia na parede deste último furo. 


c) Colocação do varão. No furo de diâmetro d é 
introduzido o varão destinado a reforçar o 
maciço rochoso (Fig. 1c). O varão é colocado 
por meio de hastes, convenientemente guiadas, 
que permitem definir-lhe a orientação. Para ligar 
o varão ao maciço é injectado aglutinante, atra- 
vés das hastes e do próprio varão, que também 
é furado (Fig. 1d). A ligação do varão às hastes 
permite que estas sejam retiradas antes de o 
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aglutinante endurecer, o que é necessário em 
virtude de o aglutinante poder envolver a peça 
de ligação e portanto dificultar a actividade das 
hastes após endurecimento. 

Têm sido usados varões furados de aço e, por 
vezes, de outros materiais, com um diâmetro 
diferindo de alguns milímetros de diâmetro do 
furo, por exemplo 5 mm. Como aglutinante têm 
sido usadas sobretudo caldas de cimento mas 
também foram experimentadas resinas epóxido. 
A fim de melhorar a sua resistência, utilizam-se 
caldas de cimento com uma relação água-cimento 
tão baixa quanto seja compatível com o escoa- 
mento do aglutinante no espaço entre o varão 
e a superfície do furo. É corrente usar relações 
de cerca de 0,4. O aglutinante é colocado num 
reservatório (Fig. 2) donde sai, através de tubo 
de plástico ligado às hastes de colocação do 
varão, mediante a aplicação de ar sob pressão. 
Quando se deseje evitar que o aglutinante pene- 
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tre no maciço rochoso é aconselhável ir levan- 
tando o varão à medida que o furo vai sendo 
preenchido e depois levar o varão lentamente 
à sua posição. Convém usar aditivos para regu- 
lação do tempo de endurecimento do agluti- 
nante, com vista à economia da operação de amos- 
tragem. Também é por vezes conveniente adicio- 
nar corantes a fim de ser facilitada a distinção 
entre o aglutinante e o restante material quando 
da observação dos tarolos. 

A fim de evitar que nos espaços preenchidos 
com o aglutinante entre água do terreno que 
prejudique a ligação do varão, é aconselhável 
que no decorrer do endurecimento a altura da 
água no furo de sondagem não seja inferior 
à correspondente à pressão da água intersti- 
cial do maciço rochoso na zona da amostra- 


gem. 


d) Sobrecarotagem (Fig. 1 e). Após o aglutinante 
ter atingido o grau de endurecimento necessário 
procede-se à sobrecarotagem da zona armada 
do maciço rochoso. É conveniente organizar as 
operações de modo que o endurecimento se pro- 
cesse no decorrer da interrupção nocturna dos 
trabalhos de sondagem. A sobrecarotagem deve 
ser feita com os cuidados habituais, sendo con- 
veniente usar coroa de pequena espessura a fim 
de ser maior a espessura da amostra. Têm 
sido obtidos resultados inteiramente satisfatórios 
mesmo usando extractor simples. No caso de 
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maciços de muito má qualidade pode ser acon- 
selhável a utilização de amostrador duplo. 
Quanto aos reforços a que se encontra subme- 
tido o aglutinante que preenche o espaço entre 
o varão e o maciço rochoso, nota-se que são em 
regra diminutos pois as acções que o extractor 
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Fig. 3 


exerce sobre o maciço se concentram num bolbo 
tal como o indicado na fig. 3. O aglutinante terá 
sobretudo de suportar solicitações secundárias, 
em especial as provenientes de vibrações, a não 
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ser quando ocorrem fendas no maciço rochoso 
com certas orientações. Assim, no caso de fenda 
aberta perpendicular ao varão (Fig. 4) as solicita- 
ções aplicadas pela coroa podem transmitir-se 
em grande parte à ligação do varão ao maciço 
rochoso. Além disso, quando a coroa se aproxima 
da fenda pode destacar-se, como já temos obser- 
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vado, o volume de rocha indicado na fig. 4, em 
consequência do desenvolvimento de tensões de 
tracção. Para atenuar este fenómeno poderá ser 
necessário reduzir a força axial aplicada na 
sobrecarotagem. 


3. APLICAÇÃO DO MÉTODO EM LABORA- 
TÓRIO 


O método foi aplicado em laboratório com 
vista a julgar a sua viabilidade e a definir as 
melhores condições da sua utilização. 

Para isso foram simulados maciços rochosos 
mediante a sobreposição de camadas de rochas 
de diversas naturezas e com espessuras de 1 a 
20 cm. Entre essas camadas eram deixados inter- 
valos com dadas espessuras, variando desde O a 
10 mm, os quais em certos casos foram deixados 
abertos e noutros preenchidos com materiais 
siltosos e argilosos de diversas consistências, 
reproduzindo os enchimentos que ocorrem na 
Natureza. Essas camadas eram montadas em 
tanques (Fig. 5) e depois procedia-se à extracção 
de amostras integrais, mediante a abertura de 
furo para alojamento do varão, a introdução 
deste e do aglutinante, e a sobrecarotagem. As 
amostras eram colhidas perpendicularmente às 
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camadas, por ser esta a direcção mais desfavorável. 

No respeitante à operação de amostragem 
foram numerosos os aspectos estudados : tipo de 
coroa e extractor adequado à abertura do furo 
para o varão, material do varão, intervalo entre 
o varão e o furo, centragem do varão, tipo de 
aglutinante, em particular a sua consistência, 
velocidade de introdução do varão no furo 
preenchido com aglutinante, aditivos para regu- 
lação da evolução do endurecimento e para colo- 
ração do aglutinante, e finalmente o tipo de coroa 
e extractor para a sobrecarotagem e condição da 
sua utilização. A técnica que aconselhamos para 
a aplicação do método, exposta na secção 2, é o 
resultado destes estudos. 

Na fig. 6 apresenta-se, na parte superior, 
amostra integral obtida em laboratório, em furo 
com diâmetro 7,5 cm aberto com extraçtor simples. 
A solidarização foi obtida por meio de varão de 
plástico e calda de cimento. Podem observar-se 
as camadas de rochas de diversa natureza e o 
seu diferente espaçamento. Os intervalos estavam 
preenchidos com materiais argilosos, verifican- 
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do-se que mesmo camadas praticamente isentas 
de coesão se mantiveram em grande parte na 
amostra integral. Para comparação, apresentam-se 
na parte média da fig. 6 os tarolos obtidos num 
furo com o mesmo diâmetro aberto junto do 
anterior, com a técnica corrente de amostragem, 
com extractor duplo. A percentagem de recupe- 
ração, determinada como é habitual, foi de 89º/,, 
não se tendo conseguido praticamente recuperar 
nenhum dos enchimentos e, além disso, algumas 
das camadas de rocha sofreram, como se vê na 
figura, fracturação apreciável. Na parte inferior 
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da fig. 6 vêem-se ainda os tarolos obtidos na 
abertura do furo, côm diâmetro de 2,2 cm, para 
alojamento do varão, a qual foi feita com extrac- 
tor simples. A recuperação foi somente de 
75 “Ja. 


superfícies cortadas permitem uma «leitura » 
perfeita de todas as fracturas e seus enchimentos. 
No caso em consideração as juntas encontra- 
vam-se abertas e por isso foram preenchidas 
com o aglutinante., 
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Na fig. 7 vê-se tarolo obtido também em 
laboratório, em furo com 7,5 cm de diâmetro, 
cortado longitudinalmente, por meio de disco 
diamantado, ao longo de duas secções radiais, 
como se está fazendo habitualmente para melhorar 
as condições de observação. Na verdade, as duas 
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4. APLICAÇÕES NO CAMPO 


Após o estudo do novo método em laboratório 
foi empreendido um programa de aplicações no 
campo, sobretudo em maciços rochosos que 
exibem condições mais desfavoráveis, com vista 
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a definir as melhores técnicas de aplicação tendo 
em atenção a economia da amostragem. 

O método foi aplicado em furos com diâmetros 
desde 7,5 cm até 17 cm, em regra em furos 
verticais, já tendo porém sido usado em furos 
com a inclinação de 45º, e quanto a profundidade, 
apesar de ainda só termos tido oportunidade de 
colher amostras até 15 m, a experiência alcançada 


condições tais que o método comum de amostra- 
gem conduzia a percentagens de recuperação 
muito baixas ou não permitia mesmo obter tarolos. 

As figs. 8 a 11 referem-se a formações grani- 
ticas decompostas ou mesmo já convertidas em 
solo residual. Apesar do estado de alteração foi 
possível obter sempre amostras de todo o maciço 
rochos». No caso da fig. 10 tratava-se de zona 
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mostra que as técnicas desenvolvidas podem ser 
aplicadas a qualquer profundidade. Para aprecia- 
ção dos méritos do novo método ele foi várias 
vezes aplicado em furos a pequenas distâncias, 
da ordem de 0,5 m, de furos nos quais a 
amostragem tinha sido realizada seguindo a técnica 
comum, com extractor duplo. 

Apresentam-se a seguir referências sumárias 
a amostragens feitas em maciços rochosos com 
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de granito decomposto onde se cruzavam três 
diaclases com enchimento argiloso muito mole 
que não permitiu a entrada do aglutinante nas 
diaclases. Apesar de condição tão adversa, só se 
perderam algumas pequenas partes superficiais 
da amostra. 

Nas figs. 12 a 16 apresentam-se alguns aspec- 
tos de amostras integrais extraídas de maciço 
constituído por um xisto clorítico, com graus 
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de alteração variáveis, que continha frequentes 
intercalações lenticulares de mármore e também 
cavidades resultantes da dissolução dessas inter- 
calações. Além disso, o maciço era cortado por 
várias falhas. Nas figs. 12 e 13 vêem-se bem as 
cavidades provenientes da referida dissolução e 
compreende-se que a percentagem de recuperação 
pela técnica comum fosse baixa. 


5. CONCLUSÃO 


Passou-se pois a dispor de método de amostra- 
gem que permite a observação de todas as carac- 
terísticas dos maciços rochosos e da sua orienta- 
ção. Mais do que isso, as amostras colhidas 
podem ser ensaiadas para determinação das pro- 
priedades das fracturas de vários tipos, incluídos 


Fig. 13 


Fig. 14 


As figs. 14 a 16 mostram como o método da 
amostragem integral permite mesmo a perfeita 
amostragem de enchimentos de falhas, sendo de 
destacar o caso da fig. 16 que se refere a amos- 
tra que incluía troço com cerca de 20 cm, de 
um enchimento argiloso mole. 
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os seus enchimentos, e das zonas de fraca resis- 
tência, em especial das propriedades mecânicas. 
Deste modo, a amostragem integral contribui 
para a aplicabilidade das teorias da mecânica das 
rochas que consideram o papel relevante desem- 
penhado pelas fracturas, pois até agora se tem lu- 
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tado com a dificuldade de as detectar e caracte- 
rizar. 

É oportuno referir que se, por um lado, a 
amostragem comum fornece informação insufi- 
ciente e conduzindo por vezes a apreciação opti- 
mista em virtude de não estarem representadas 
as zonas de menor resistência, por outro lado, a 
observação das paredes de galerias de inspecção 
pode conduzir a apreciação desfavorável, dada a 
perturbação do maciço rochoso derivada da aber- 
turá das galerias. 

Nota-se que, enquanto na técnica comum de 
amostragem é em regra difícil verificar que o 
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grau de recuperação obtido é o que séria razoável 
exigir, ma amostragem integral a situação sim- 
plifica-se porque só é aceitável, salvo casos de 
espécie, uma recuperação de 1009/,. 
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NOTICIÁRIO 


Preparação do 5-º Congresso do Conseil 
International du Bátiment 


—e 


A preparação do 5.º Congresso do CIB, que 
terá lugar no Palácio dos Congressos de Ver- 
sailles de 22 a 30 de Junho de 1971 prossegue 
activamente. 

A orientação geral do Congresso, cujo título é: 
«Da investigação à prática: o desafio da aplica- 
ção», é a de constatar o estado dos conheci- 
mentos úteis à construção e o grau da sua apli- 
cação na prática. Espera-se que será a ocasião 
para as pessoas ligadas à construção de melhor 
conhecer o que já existe e também para deter- 
minar o que seria necessário fazer para que o 
que se conhece seja mais divulgado e utilizado. 

Praticamente todo o domínio das ciências e 
técnicas da construção será coberto pelos 16 
temas do Congresso : 


1 — Os dados naturais 
2 — As exigências dos utilizadores 
3 — Acústica 
4 — Higrotermia 
5 — Informática 
6 — Estática 
7—0O fogo 
8 — Economia da construção 
9 — As ajudas para a concepção 
10 — Tecnologia do betão vazado 
11 — Sistemas de construção ligeira 
12— A tecnologia industrializada das casas 
individuais e das escolas 
13— À tecnologia dos prédios altos 
14 — Tendências no emprego da madeira 
15 — Novas tendências no equipamente e ven- 
tilação 
16 — O habitat muito económico. 


Uma brochura dando todos os dados sobre os 
temas, o calendário, o local das sessões e as 
informações práticas de toda a ordem pode ser 
obtida no 


Sécretariat du Congrês du CIB 4, Av. du 
Recteur Poincaré 
Paris, 16 eme. 
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COPRAI 
Aperfeiçoamento dos métodos e estudo 
dos tempos na Indústria de Artes Gráficas 


A COPRAI vai organizar, de 29 a 31 de Março 
de 1971, um curso sobre o tema indicado com a 
duração de 18 horas. 

O Curso que será ministrado em Lingua 
Espanhola pelo Eng. D. ENRIQUE DEL NOGAL 
ROPERO terá o seguinte 

PROGRAMA: 

Produtividade 

Estudo do trabalho 

Produção 

Aperfeiçoamento de Métodos 

Diagramas 

Implantação 

Estudo de Movimentos 

Exercícios práticos, Investimentos e aperfeiçoa- 

mentos tecnológicos 

Estudo dos Tempos 

Estudo dos Tempos com cronómetro 

Valoração da actuação 

Suplementos 

Tempo padrão 

Simplificação de trabalho administrativo. 


Feiras e Exposições em Basileia em 1971 


No ano de 1971 realizar-se-ão em Basileia as 
seguintes feiras e exposições interessando à 
engenharia: 

— De 9 a 13 de Março 

5.º Salão Internacional de Electrónica Indus- 
trial — INEL 

— De 17 a 27 de Abril 

55.º Feira Suiça de Amostras 
— De 8 a 12 de Junho 
5.0 Salão Internacional da protecção do 
meio vital 
— De 3 a 23 de Junho 
Exposição industrial da URSS 

— De 13 a 18 de Setembro 
5.0 Salão Internacional da Técnica de Labo- 
ratório, da Técnica de Medida, de Engenha- 
ria Química e da Automática em Quimica 

— De 24 de Setembro a 2 de Outubro 

HOLZ 71 — Feira para o trabalho da Madeira 
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INVERSÃO NUMÉRICA DE TRANSFORMADAS DE LAPLACE 


RESUMO 


Sugere-se um processo numérico para inversão de um 
tipo restrito de transformadas de Loplace e refere-se a 
sua aplicação imediata no domínio das técnicas de 
regulação. 


1. INTRODUÇÃO 


As transformadas integrais desempenham, nos 
problemas técnicos, um papel de crescente impor- 
tância, permitindo a manipulação de questões 
insolúveis por outras técnicas. A transformada 
de Laplace reúne um potencial notável, sob esse 
aspecto, sendo até a mais amplamente divulgada. 

A transformada de uma função f (t) define-se 


dp=“(t0)=[ te" a 
“o 
e a sua inversa 


E(t) = 47 E (p)| na 


2+joo 
= [ ó(p) e” dp 


E — | DO 


A determinação analítica da transformada 
inversa envolve, por vezes, operações delicadas 
ou, pelo menos, laboriosas. Estas poderão ser, 
nalguns casos, torneadas, utilizando processos 
numéricos de cálculo. No presente artigo propõe-se 
um método numérico abrangendo, apenas, um 
sector restrito das transformadas inversas. 


2. TRATAMENTO MATEMÁTICO 
2.1 — Posição do problema 


Procura-se a determinação de transformadas 
inversas do tipo 


É (t) = Lo 2, 
Pp 


com p=«-iw, respeitando-se as condições 
seguintes : 


1) 2==0. Será então 


F(p) = F(io) = Ulw) +i V(v) 
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À numerical method for inversion of a restricted 
range of Laplace transforms is sugested and its imme- 
diate interest in the field of regulation technics is referred, 


2) F(p) é uma função racional de coeficientes 
reais, tal que 


F (— e) = U (1) ais V (mm) 
2.2 — Tratamento formal 
O leitor interessado poderá encontrá-lo em [1]. 
Aqui transcrevem-se, apenas, as expressões que 


permitirão solucionar o problema enunciado 
em 2.1. 


sen wt dv 


É (t) = 


ar 


2 f <a 


o 


E(t) = + [ V (m) cos mt do + Ele) 
Tr “ fa 
0 


As duas expressões podem ser utilizadas inde- 
pendentemente uma da outra e ao critério do 
utilizador, apenas havendo a notar que: 


U (1) V (13) E : 
1) ———— ou ——— terão que ser tais que 
ts tm 


lim + ta) = 0 lim o tes. () 


Q — CO ty 


OD > Cx fm 


2) Para uma mesma precisão de resultados, 
exigirá menos tempo de cálculo a integranda 
mais rápidamente convergente para zero, o que 
se pode avaliar através dos 


o U(») 
lim —— 
us co V (1) 

3) Para se utilizar U »', terá que ser U o) 
finito ou nulo; para utilizar V(»m) deverá ser 
V(0)=0. 

4) U (») ou V(») não terão zeros para valores 


de vw =0. 


As condições impostas em 21—1) e 2) e 
2.2 — 1), 2), 3) e 4) não restringem drasticamente 
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o campo das funções a que o método é aplicável, 
posto que muitas são as que as satisfazem. 


2.3 — Tratamento númérico 


Para cada valor de t será possível representar 
gráficamente qualquer das integrandas. Vamos 
supor a representação fictícia da Fig. 1. Os zeros 
da integranda serão os zeros da função sen »t. 
É evidente que a integração não poderá ser 
levada até ao infinito e deverá, portanto, ser 
interrompida em certo »,,. É também óbvio que 
o erro cometido será menor se »,, se situar entre 
dois zeros e será ideal para certo valor entre os 
dois zeros escolhidos Não se procurará o valor 
ideal, mas sim o mais simples e esse será o ponto 
médio entre dois zeros consecutivos, facílimo de 
determinar em funções do tipo sinusoidal, Os 
valores absolutos da função integranda nos pon- 
tos médios serão, sucessivamente, comparados 
com uma quantidade tão pequena quanto quiser- 
mos, satisfazendo o que tiver valor menor ou 
igual a essa quantidade. Na figura 1 as linhas 
ponteadas representam o valor absoluto da inte- 
granda nas zonas onde ela é negativa e a linha 
horizontal o valor, por assim dizer, padrão. No 
caso da figura, », seria 2 3.º ponto médio e a 
integração processar-se-ia de 0a pela regra de 
Simpson. Nesta, o erro em seádão ao valor real 
do integral será controlado pela diferença — ER- 
RO2 — entre dois integrais aproximados suces- 
sivos achados com um número de intervalos 
duplo um do outro /2]. A precisão, controlável 
pelo valor ERRO 1 e ERRO 2 será tanto maior 
quanto menores eles forem, mas o tempo de cál- 
culo crescerá, também, simultâneamente. O racio- 
cínio, para o caso de se usar V(»), é perfeita- 
memte paralelo ao exposto. 


2 4! 4 
Km senqw- 
Ea Uli) sen w 


ta, 
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2.4 — Exempios 


24.1 — Considere-se a função 


1 
E(p) = —— 
1+p 
que poderá escrever-se 
| fa) 
É 5) == ———— ms. comi i - 
nd pr Da E 
U (si) 


O valor de lim permite concluir ser 


to —s OS Es 


U (») mais favorável para os nossos propósitos. 
No entanto, e porque se trata de um exemplo de 
aplicação do método, utilizaremos, paralelamente, 
U(m») e V(») para o cálculo de f(t). Reprodu- 
zem-se tabelas (TAB I e TAB II) em que se com- 
param os valores obtidos numericamente e os 
valores correctos conseguidos por cálculo exacto 
a partir da expressão f(t) == 1 — e”!, deduzida 
analiticamente por inversão 


ge] [=| 
pt +p) 


Nas figuras 2 e 3 reproduzem-se, também, os 
programas que permitiram realizar, num compu- 
tador IBM 1620, os cálculos descritos. Os pro- 
gramas poderão ser usados para qualquer fun- 


POQ FOSTRAN C2 

START 

- ES bOS UIMENSION tada FOM(15),5415] 

- GÊ U | ça rir 28, 1) 

532 é CAMATÇF 7,4 li 

-D Sd - ECRMAT IS] 

-sÊ 76 SU FORMATIF7,3,2E14,8) 
-6706 READ &,N 

3734 REÃO 1,PI,ERRDOI ERRO? 
-t Tel DO 50 1=1,N 
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ção, adaptando ao caso em estudo os «state- 
ments» assinalados com uma seta, que represen- 


tam, por ordem, de cima para baixo 


a —bT)o—b,T, wo 


RT at) +41. 


com 
16000100000M 

START 

-6600 DIMENSION T(15),FOM(15),S(15) O valor de 

-—6600 1 FORMAT(F11.6,2E6.1) 

-6632 2 FORMATIET 3) : 

-b654 by FORMAT 3) lim UÚ (19) 

- N f 

673? READ dd sor Eras % —s OO o) 
6792 50 READ 2T(Ã) À ? 
6976 DO 400"1=1,N mostra ser mais conveniente para os nossos 
“6900 70 OM=V*PI/T(I) cálculos U(»w). O novo programa, com os 
“6972 FOM(1)=-2.*C0S(OM*T(I))/(PIR(I + OM AZ) | cr a o a à 
7164 rd rd PRA «statements» já assinalados adaptados à nova 
-77280 LO Vev+!. E a 

[5516 GO TO" 70 função, reproduz-se na fig. 4. Numa tabela 
++ Aa (TAB III) comparam-se os valores obtidos 
-7372 FO=-2./P| eim Es 

-J42 SO=FOIFOM( 1) numeéricamente com os exactos, calculados a 
= TU! = 

-7492 Sil==2, *COS(OX*T(I))/(PIN(1. COXA) |) “ “ a fas 

EE SUTIL ASONR AS IDR. “po partir da inversão analítica 

77 104 ( 

7792 LegAl 

-7828 Kel =] Pas b, + b, PU 8 A 

- 186º bio & 

— Õ = “5 | 

29935 S 100, P(Lp pP+HTp+1) 

-19h8 DO 2C0 Jet,K,2 

25398 209 to S 1-2 .*C0S((2*DK)*T( Ci (Z*DX)w*2)) ar 

- 14138 ERRA = — *Cco * + + * -—— —= 

-8736 S(1)=( Si . 85 102 .45P )4DX 3. CONCLUSÕES 

- 3300 A=ABS(S(1)-S1) 

-B452 IF(A- ço 306,100,100 


-8520 300 5(1)=S(1)+ 1) O processo numérico couduz a valores per- 
- SO; 400 PRINT 60, Ti), FOMCI),S(I) 


Rr ENO. feitamente aceitáveis para as necessidades de 
precisão correntes, podendo tornar-se mais satis- 


PROCESSING COMPLETE fatórios ainda à custa de um aperto dos valores 


START . . 
ERRO 1 e ERRO 2, implicando, todavia, um 
Fi * e 
8. 3 PDO FORTRAN C2 
START 
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rw T (1 -6676 80 FORMAT(F7,3,2714.8) 
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para ww == OM -6732 READ 1,PI,ERRDI ERROZ 
Id eae 50 REMO 2 Tli) 
2) Idem para »=0 6876 DO 300 I=1,N 
as & V=1. 
3) Idem para w = DX 6900 70 Ona(zemyss asi 42 SI) 
, E =2.*(12.5+262 .6%0M4%2 )*SIN(OM*T( 
4) Somatório dos termos ímpares em que ig ETAA VEMOS AS BLOOD ADIAR 
2U (0) =1332 FOM(1)=B/C 
a 2.2 parcela será o valor da função — Elisa pr ibid 0,10,40 
; 6) -7508 40 VeV+l. 
- 7 Sul GO TO 70 
t 2 V (0) Z> DX Ju O oco. sum 
(5) ———— (1) - X=0. 
sen ou mn cos wt| para >< “3800 FO=2.*12.5*T(1)/PI 
E) -76Bh S0=FO+FOM( 1) 
cê Bi=2 al t2 55262. GeDNMAD) AS INC ORAT( 1) 
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Notar que, no caso de se utilizar V(w), o «sta- 78! 2 Cla(ão. 322Dxtebiiat. B6*DXART 4, gut 2x 
. . ss = 1 
tement» 300 adiciona ao integral aproximado  -505: (1 )=(S0+b 451 )4DX/3. 
S(1) o valor F(0) Ea 
de Le2* 
Valor -2260 K=L-1 
-2296 DX=0.5*0X 
24.2 — Tome-se a expressão -8342 SP=SPSi 
8350 DO 200 Ju1,K,2 
b.i+b (io) -8302 Z=) 
E (i (à) o o l 31428 B2=2,*(12,5+252, io ger Jo À AS 3) 
ENA E T (i 0)? E T (i w) + 1 Ee] E C2=(40.324(2%0X) tt 121, 86*(2*0X)**2+7.)nPI*X(2*DX) 
- = EE 
2 l -8896 S(1)a gpa 45 142 .*SP)*DX/3, 
-30h40 A=ABS(S(1)-S1) 
ci -9112 IF(A-ERROZ) lo, 100,100 
-9180 300 PRINT 80,T(1),FOM(1),S(I) 
2 25358 END. 
bb+(b T;—b, T;) PROCESS ING COMPLETE 
+ START 


F (00) um 


*+(T-2T) o +1 Fig. 4 
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inevitável aumento do tempo de cálculo. Não há 
um critério absoluto para a escolha destes dois 
valores: apenas a experiencia, a análise do tipo 
da função em causa e a precisão pretendida per- 
mitirão a selecção mais adequada (!). 

2) O processo tem aplicação imediata nas 
técnicas de regulação envolvendo sistemas lineares 
em que, como se sabe, a resposta harmónica se 


define 


f Lf 
E (p)= LF) ss Lito, 
L(s(t)) 1/p 


em que s(t) será tal que 


s(t)=0 tão 
s(t)=0 tãú 
e p= im». Então 
| Fíp) 
E (1) =| FE 
LP 


3) O método poderá prestar algum auxílio em 
problemas envolvendo transformadas de funções 
satisfazendo aos requisitos enunciados. 


TABELA 1 


Valores de f(t) usando 


U (») = —, com ERRO 1=ERRO 2=0,001 
44 + tw) 
E(t) E(t) 
Processo numérico | Cálculo exacto 
| 
0,449 | (1.3622 0.3617 
0,628 | 0,4062 0,4663 
0,785 | (1,5436 (1,5439 
1,047 | 0,0491 0,0490 
1,570 1,7018 11,7919 
3,140 (1,95307 0,U567 
3.950 0.08U8 O 0807 
| 6,280 0,098] 1,998] 


() Em casos correntes, considera-se aceitável o valor 
ERRO 1 = 107º, podendo ser aumentado quando o poli- 
nómio denominador for de grau elevado. Para ERRO 2, 
um valor da ordem 107? garantirá aproximações normal- 
mente satisfatórias. Os valores indicados devem ser 
tomados como base e não como universais. 
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TABELA II 


Valores de f(t) usando 


4) 


V(m) =— 5 com ERRO 1=ERRO 2=-0,001 
1 = fa) 
| E(t) E(t) 
' Processo numérico | Cálculo exacto 
| | 
(1,449 (1,3620 0,3017 
0,028 | 0,4065 0,4003 
(0,785 (1,5438 | 0,5439 
1,047 0,0489 (0,0490 
LáTo 0,7010 O,7919 
3.140 | 0,9567  0,9567 
3,950 09807 | 0,9807 
| 6,280 | 0,998] | 0,998] 
| | 
TABELA III 
Valores de f(t) usando 
bo + (b,T, - bo T,)0º 
AE ES A e SD , com 
Tin +(T;—-2T)w+1 
ERRO 1 = ERRO 2=-0,01 
E (t) | E (t) 
' Processo numérico | Cálculo exacto 
0,449 Lo ame 1,562 
(0,785 2 2080 | 23005 
VSTO 3,3007 33705 
340 4,0946 | 4.044 
(0,280 | (1,0354 | (,0308) 
8,30) 1055 FTOSO 
12.500 | 0,1045 | 01646 
16,500 10,1795 10,1795 
28 000 11,0824 11 6824 
45,000 | 12,325] 12 325] 
| 60,000) | 12,455] 12,455] | 
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Pontes rolantes 
fabricadas em cadeia... 


As pontes rolantes Standard são fabricadas em série, nos 

tipos monoviga e biviga, pelo sistema de combinação de 

peças normalizadas. 

Para vãos até 28m e capacidades de carga de 250 kg 

a 20.000 kg permitem a adaptação às vossas condições 

de trabalho. : 
Reduzido peso próprio, pequena altura de construção, 
prazos de entrega rápidos, preços excepcionalmente 
baixos, eis as vantagens das pontes rolantes standard 
«Demag»: 
Consulte-nos e os nossos engenheiros efectuarão o es- 
tudo do seu problema gratuitamente, 
Assistência exemplar - O Vosso sucesso será também o 
nosso. 
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DESCRIÇÃO DE SISTEMAS NOS DOMÍNIOS DO TEMPO, 


IMAGEM E COMPLEXO 
DE TRANSFERÊNCIA E 


RESUMO 


Define-se o conceito e tipos de sistemas e efectua-se 
a respectiva descrição através da equação diferencial, 
da função de transferência complexa e da resposta de 
frequência. A partir das equações dinimicas define-se 
a matriz de transferência complexa e caracterizam-se 
os sistemas por meio da equação de estado, mostrando-se 
as propriedades mais pertinentes e menos csclarecidas 
na literatura, 


1. INTRODUÇÃO 


Na técnica de controle utiliza-se a noção de 
transferência para se caracterizar a resposta a 
uma determinada excitação, em geral através 
da função de transferência complexa do sistema. 
A aplicação do método é de tal modo frequente 
que se perde ocasionalmente a medida das suas 
limitações, o que provoca uma certa imprecisão 
nos estudos por vezes apresentados na literatura. 

Por exemplo, o estabelecimento da equação 
diferencial de um sistema não é normalmente 
caracterizado, conduzindo a uma imperfeita, às 
vezes até falsa, noção de função de transferência 
complexa. Igualmente se encontram misturados, 
chegando-se a afirmar a sua identidade, os domíi- 
nios imagem e complexos, 

Pretende-se neste estudo caracterizar as noções 
essenciais da teoria dos sistemas, fundamento 
da teoria da automação, por meio dos seus 
conceitos básicos e propriedades fundamentais. 


2. CARACTERIZAÇÃO DE SISTEMAS 


Chama-se sistema a todo o conjunto de elemen- 
tos físicos, caracterizado por um determinado 
número de parâmetros, que pode transferir 
qualquer sinal de entrada (aplicado aos terminais 
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ATRAVÉS DO REGIME 
DO ESTADO DINÂMICO 


HERMÍNIO DUARTE-RAMOS 


Engenheiro Electrotécnico (1.S.T.) 
Colaborador Científico do Shering-los- 
titut (T. U. Hanover) 

Bolseiro da U, L. (Angola) 


por 


SYNOPSIS 


The concéption of system and its types are defined 
and particular systems are described by differential equa- 
tion, transfer function and frequency response. The dy- 
namic equations and the transfer matrix are introduced. 
The state equation is used to define the principal sys- 
tems, and the most important properties are criticised, 


de entrada do sistema) para os terminais de saída 
(onde se capta o sinal de saída), ou, mais genê- 
ricamente, é qualquer estrutura física que permite 
fornecer uma informação. 

Os parâmetros de um sistema são grandezas 
físicas que caracterizam os elementos físicos de 
constituição do sistema, por exemplo, a capaci- 
dade de um condensador ou a massa de um corpo. 

Se o valor dos parâmetros de um sistema 
depende da posição no espaço e do tempo (é o 
caso de uma linha de transmissão) diz-se que 
se trata de um sistema de parâmetros distribuídos ; 
se os parâmetros são apenas função do tempo 
(e não do espaço) trata-se de um sistema de pará- 
metros concentrados. Em geral a aproximação dos 
sistemas de parâmetros concentrados é suficiente 
tornando-se mais simples a sua análise. 

Na definição de sistema refere-se o conceito 
de sinal que é um estado físico utilizado para 
reconhecimento de uma informação (sinal de 
entrada ou excitação do sistema) e doação dessa 
informação (sinal de saída ou resposta). O sinal 
de entrada equivale a uma excitação ou causa e 
o sinal de saída a uma resposta ou efeito. Os 
sinais, sendo observáveis apenas através de varia- 
ções no tempo, são funções do tempo, represen- 
tando-se, convencionalmente, o sinal de entrada 
por v(t) e o sinal de saída por x (t). 
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Esquematicamente pode representar-se um sis- 
tema por meio de um esquema de blocos onde figu- 
ram os sinais de entrada e de saída (fig. 1). 


Causa Efeito 
yce Excitação Resposta xe) 
d) 
Letts | R ct) 
uct QUADRIPOLO E vact) 
b) 


Fig. 1 — Esquema de blocos de um sistema. 
a) Caso geral; 
b) Quadrípolo. 


Um sistema é em geral constituído por certo 
número de elementos concentrados armazenado- 
res de energia, como sejam condensadores e 
solenóides em circuitos eléctricos ou molas e 
e massas em circuitos mecânicos, dizendo-se 
portanto que se está em presença de um sistema 
com memória. 

O momento em que se começa a observar o 
comportamento de um sistema denomina-se ins- 
tante iniçial de observação t,. Este valor t, é 
arbitrário e fixo e corresponde normalmente ao 
momento em que se aplica o sinal de entrada 
(mas pode ser outro valor, como é forçoso no 
caso em que não existe sinal de entrada). 

Se os elementos armazenadores de energia se 
encontram completamente descarregados no ins- 
tante inicial (para t==t,) então refere-se a um 
sistema sem energia. Este é o caso mais importante 
na prática, devendo contudo observar-se a exis- 
tência das suas condições a fim de evitar falsos 
resultados. 

Um sistema pode não ter qualquer exci- 
tação à entrada, caracterizando-se assim um 
sistema livre 


y(t)=o0 parat >t, (t, qualquer e fixo) (1) 


onde t, é o instante inicial de observação. 

O princípio de causalidade pode agora enunciar-se : 
num sistema sem energia não há efeito (sinal de 
saída) sem causa (sinal de entrada) 
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x(t)=0 para y(t) O et.>t, (t,qualquerefixo(2) 


Deste modo sistema causal é aquele em que as 
grandezas de saída (após o instante inicial: para 
t > t,) dependem apenas da história do sistema 
(anterior ao instante inicial: para t t,) e das 
grandezas de entrada (passado e presente), tendo 
o futuro o peso zero. 

Um sistema linear é todo o sistema onde se veri- 
fica a lei de sobreposição, correspondente à sobre- 
posição física de sinais, 

; n n 
(D50)=DS to 0) 


ES [==] 


e a lei de homogeneidade, fisicamente equivalente 
à amplificação de sinais, 


f (a.x(t) =a.f(x(t)) , com a = const. (4) 


ou seja, para qualquer sinal de entrada y (t) em 
qualquer instante t verifica-se 


) x; (t) = 3 a; L (y;(b] = 


i=1 | = 


| 
t4= 


» Lfa,y;(t) (5) 


sendo a, constantes e L [...l um operador ou uma 
transformação. 

Sistema de tempo invariante, também designado 
sistema estacionário, é aquele cujo comporta- 
mento é independente do instante de observação 


x(t—2=Liy(t—u) , 


para todo t e Z qualquer e fixo (6) 


Chama-se sistema autónomo a um sistema de 
tempo invariante que seja livre (não tendo por 
conseguinte sinal aplicado à entrada). 

Todo o sistema causal em que para grandezas 
reais também as grandezas de saída são reais 
diz-se sistema realizável. Um tal sistema exige o 
conhecimento dos sinais de entrada y (t) e de 
saída x(t) no intervalo (— co,t) ; quando apenas 
se conhece o sinal de entrada no intervalo (t,,t) 
é necessário o conhecimento das condições iniciais 
para a sua completa descrição. 

As condições iniciais de um sistema são os valores 
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que caracterizam a descrição do sistema no ins- 
tante inicial de observação t,. Estas condições 
iniciais são função dos valores dos parâmetros 
do sistema e da interligação dos correspondentes 
elementos físicos no instante inicial. Como é 
evidente as condições iniciais de um sistema sem 
energia são todas nulas. 

Em relação à duração do sinal de saída dis- 
tinguem-se : sistema contínuo (fig. 2-a) se a saída se 
processa continuamente no tempo, sistema discreto 
(fig. 2-b) no caso da saída ser definida sômente 
em certos valores discretos do tempo e sistema 
discontínuo (fig. 2-c) se a resposta se manisfesta 
por intervalos de tempo. 


Cc) 


Fig. 2— Sistemas de trabalho contínuo (a), dis- 
creto (b) e discontínuo (c). 


Todo o sistema determinista é representável por 
uma função do tempo, tendo portanto tratamento 
matemático, enquanto sistemas não deterministas, 
relativos a sinais não controláveis, apenas poderão 
ter um tratamento estatístico. 

Neste estudo considerar-se-ão somente siste- 
mas de parâmetros concentrados, causais e de 
trabalho contínuo. Normalmente refere-se a sis- 
temas de tempo invariante e em especial a sistemas 
lineares. 


3. DESCRIÇÃO DE SISTEMAS ATRAVÉS 
DO PROCESSO DE TRANSFERÊNCIA 


3.1. Descrição no domínio do tempo t 


O regime estacionário de um sistema, atingido 
uma vez decorrido um tempo teóricamente infinito 
após a aplicação do sinal de entrada y(t), 
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quando reinam dados valores de perturbação z (t) 
do sistema (fig. 3-a), determina um certo conjunto 
de valores de saída x(t). Fisicamente o regime 
estacionário é atingido decorrido um intervalo de 
tempo finito considerado suficientemente longo 
para que ao valor de entrada corresponda de 
facto o respectivo valor estacionário de saída. 


Lz(t) 


b) C) 


Fig. 3 — Regime estacionário de um sistema. 
(a) Esquema de blocos; 
(b) Superfície característica; 
(c) Família de características. 


O conjunto de pontos de funcionamento esta- 
cionário define a superfície característica do sistema 


(fg. 3-b). 
(7-a) 


(7-b) 


Y (x (t), y (t), z(t) =0 
ou x(t)=g (y(t), z(t)). 


Os valores particulares x,, Yo, Z; determina- 
dos por um dado funcionamento do sistema em 
regime estacionário caracterizam o chamado 
ponto de trabalho do sistema, 

Como a superfície característica tem em geral 
uma forma qualquer e o ponto de trabalho pode 
ser qualquer é importante o conceito de lineari- 
zação (se possível), que se obtém pelo desenvol- 
vimento de Taylor da superfície característica 


em torno do ponto de trabalho 
dg dB), 
Ax=([-2|] A ms Á 8 
x ( ) dy+ (52) dz (8 


com os desvios àx = x — x, , dy =y—y, 
e Sáz=7— 2. 
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Normalmente toma-se z(t) como parâmetro 
reduzindo-se a superfície característica a curvas 
no plano y(t), x(t). Para cada valor particular 
da grandeza de perturbação z (t) estabelece-se 
então a chamada característica do sistema através 
da curva, em geral não linear, que descreve o seu 
regime estacionário, isto é, através da relação da 


resposta x (t) a uma dada excitação y (t) quando 
t mto (O) 


x (t) = (y (t)) , com z (t) = const. (9) 


Esta função é em geral apenas representável 
graficamente. 

Ao conjunto de curvas características relati- 
vas ao funcionamento estacionário do sistema 
para distintos valores de perturbação chama-se 
família ou campo de características (fig. 3-c). 

A linearização consiste portanto em substituir 
na vizinhança do ponto de trabalho a curva ca- 
racterística pela recta tangente no ponto de tra- 


balho 


Dig; = e) ày , com z (t)=z, 


(10-a) 
d Y/o 


e, efectuando uma mudança de coordenadas para 
o ponto de trabalho (x, = y, = 0) , com K = 
(dF/dy), obtém-se 


x=Ky (10-b) 


onde x e y têm no novo sistema de coordenadas 
o significado dos desvios relativamente ao ponto 
de trabalho das grandezas x(t) e y(t), isto é, 
x= 0x e y= dy. 

Genericamente o processo dinâmico de transferência 
de um sistema estabelece-se pelas grandezas de 
entrada e de saída e respectivas variações, o que 
se traduz pela equação diferencial, em geral não 
linear [1], 


(m) 


)=0m<,n (11) 


4 (nt) “ 
de 1 Y rºcosYy 


sendo portanto a linearização, em relação ao 
ponto de trabalho, 


d 
st= (0) Ax+ e Ax +... 
dx ' dx 


| / 
dpdo) AX + to iz 
0X /0 dY /o 
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e como as derivadas no ponto de trabalho são 
constantes (por exemplo dx,/dt = 0, donde 
dx '==dx/dt—dx,/dt==x') deduz-se 


e Taio ax +a,x—x)=b +... + 
+ by +bo(y—yo), (13) 


m<n 
que conduz à equação diferencial do sistema, 
equação diferencial linear de coeficientes cons- 
tantes resultante da linearização do sistema com 
uma mudança de coordenadas para o ponto de 
trabalho 


ax(D+...+a x (t) + a,x(t)=boy(t) + 


+by() +... +by",m<n 


(14-a) 


ou b> a; x (t) = ho b,. vt) ;,m<n (14-b) 


1=0 k=0 


onde n é o número de armazenadores de energia 
que participam no processo de transferência 
desde a entrada até à saída do sistema. A equação 
diferencial de um sistema esclarece apenas o 
regime dos terminais do sistema (entrada e res- 
posta) e em princípio nada diz sobre a realização 
(sobre a dinâmica interna) do processo de trans- 
ferência, quer dizer, 


a, 


E (xx... x) =R (y, Wave) (14-c) 


As condições iniciais de um sistema descrito 
pelo regime de transferência no domínio do 
tempo correspondem às condições iniciais (para 
= t,) da sua equação diferencial 
EAD ES q vg E] UAI 
(15-a) 


sendo todas nulas no caso de sistemas sem 
energia: 


x eP 1) (t,) és 4-2) (t,) ES ca (t,) : x(t,) =) 


(15-b) 


A equação diferencial homogénea corresponde ao 
sistema sem excitação à entrada, ou seja, ao sis- 
tema livre 


n 
ba ã; x (t))=0 


i=0 


(16) 


TÉCNICA N. 404 


é vasto o nosso campo de acção 


6 DELEGAÇÕES NA METRÓPOLE E ULTRAMAR 


CONSTRUCOES TÉCNICAS, LDA. 


LISNAVE — Estaleiro da Margueira — Doca n.º 10. 


OBRAS PUBLICAS 
CONSTRUÇÃO CIVIL-INDUSTRIAL E URBANA 


PROSPECÇÃO GEOTÉCNICA * FUNDAÇÕES ESPECIAIS * CONSOLIDAÇÕES 
OBRAS MARÍTIMAS E PORTUÁRIAS * PONTES E VIADUTOS * SILOS E CHAMINÉÊS 
INSTALAÇÕES INDUSTRIAIS 


TÉCNICA —- XVII 


nas é vasto o nosso campo de acção 
EE 6 DELEGAÇÕES NA METRÓPOLE E ULTRAMAR 


CONSTRUCOES TÉCNICAS, LDA. 


| 


- 


ar my? 
1 SRA 
AR ta. so de 
- ; »4 4 


T- 
4 AP 


i " 1 a , GE ça * ; % 
é e stress prt os 


MOAGEM DA BEIRA - MOBEIRA — Moçambique — Cidade da Beira — Construção de oito células 
cilindricas com diâmetro interior de 7 metros e altura de 26,5 metros (moldes deslizantes). 


OBRAS PUBLICAS 
CONSTRUÇÃO CIVIL-INDUSTRIAL E URBANA 


PROSPECÇÃO GEOTÉCNICA * FUNDAÇÕES ESPECIAIS * CONSOLIDAÇÕES 
OBRAS MARÍTIMAS E PORTUÁRIAS * PONTES E VIADUTOS * SILOS E CHAMINÉES 
INSTALAÇÕES INDUSTRIAIS 


donde resulta o polinómio característico pela substi- 
tuição das derivadas de ordem sucessiva por 
potências do parâmetro À 


QO (0) = Dad (17) 
i=0 


que determina a equação caracteristica (equação 
homogénea definida pelo polinómio característico) 


Q ()=o0 (18-a) 
ou E 

ay=o (18-b) 

i=o 


cujas raízes (valores de 4) definem os zeros da 
equação diferencial homogénea. 

Então a solução da equação diferencial homo- 
génea, para r raízes reais de multiplicidade m, e 
valores À, com s==1,..., | eparac raízes pares 
complexas de multiplicidade m, e valores 2, com 
c=pv+1,...,vserá 


+ > (exp (A, t) : > Cia t'+exp Qjtlejs!) 


s=u+l N i=1 
(19) 


sendo a solução geral da equação diferencial não 
homogénea 

X=%y + Xp (20) 
onde x, é uma solução particular, especial para 
cada caso concreto. 

Na prática o processo de transferência de um 
sistema no domínio do tempo é caracterizado 
pela resposta aos chamados sinais de teste. Nor- 
malmente usa-se a função impulsiva unitária (fig. 4-a) 
definida por 


O parat<t, 


(21) 
1 parat Dt, 


1 (t-t,) -| 


à qual corresponde a função resposta unitária, 
também designada função transição (fig. 4-b). 

Em presença de um dado sistema, caracteri- 
za-se o seu comportamento pela aplicação de um 
sinal impulsivo unitário na sua entrada e medin- 
do-se o respectivo sinal de saída. Assim se espe- 
cificam, por exemplo, os comportamentos propor- 
cional, diferencial e integral dos sistemas (fig. 5). 
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Fig. 4 — Funções de teste (a, c) e respectivas 
respostas (b, d). 
a) Função impulsiva unitária; 
b) Função transição ; 
c) Função impulso ou de Dirac; 
d) Função peso. 


x(E) x(t) 


0 t 
d4) 
x(t) 

(6) t 
b,) 
x(t) 

0 E 
C,) 


Fig. 5 — Funções transição de sistemas 
ideais (1) e com atrazo (2). 
a) Sistema proporcional; 
b) Sistema integral; 
c) Sistema diferencial, 
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Além disso interessa considerar a função impulso 
ou de Dirac (fig. 4-c) definida pelas duas condi- 
ções seguintes 


GO 


o 
fe (t) 0 (t-t,) dt = Ff (t,) 


— 


Pe 


sendo então a resposta impulso (fig. 4-d) chamada 
função peso (da teoria das distribuições). 


3.2. Descrição no domínio imagem s 


A aplicação da transformação de Laplace no 
tratamento do regime de transferência de um sis- 
tema apresenta duas vantagens de enorme valor : 
em primeiro lugar transforma a equação diferen- 
cial numa equação algébrica, mais fácil de tratar, 
e em segundo lugar o sinal de saída obtém-se 
por meio de operações apenas algébricas (utili- 
zando tabelas). Mas esta transformação implica 
a mudança do estudo do domínio do tempo t 
para o domínio imagem s=7 + j» onde, fisica- 
mente 7 tem o significado de uma constante 
de amortecimento e w de frequência angular. 

A transformação de Laplace consiste na transfor- 
mação da função f(t) no domínio t numa função 
F (s) no domínio s tal que 


F(s) = 2 (E (t) | = [oe exp(—st) dt (23) 


e a transformação inversa de Laplace, efectuando a 
passagem do domínio imagem para o temporal, 
é definida por 


E(t)=L"" TE(s)) = 
1 A 4 e oo 
Po / F(s). exp (st) ds, t >0 (24) 
cad O MR 


Nestas condições, por exemplo, a transformada 
de Laplace da função impulso de Dirac é igual 
à unidade 


“tgi ())= [05(0. exp (—st) dt=1 (25) 
como se verifica pela definição de 9 (t). 
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A equação transformada de Laplace (no domi- 
nio s) da equação diferencial de um sistema 
(definida no tempo t) chama-se equação imagem 
do sistema, assumindo a forma, para t,=0, 


a, [s"X(9) —x(+0)s"—! = x(+ 0)s"?— 


emo cem" t+ 0) | 

+a"o sro! X(s)—x (+ ve ge 0)". 
SR at |+ 

1 EP SER EE ap + (26) 


+ a,| Ss“ X()—x(+ 0)s—x'(+0) |+ 
+a,| sX(s)—x(+ 0) | + 
+a,|X(s) |= 
=|botbis+...+bys” | Y(9). 


No caso de um sistema sem energia todas as 
condições iniciais são nulas e portanto a equação 
imagem reduz-se a 


(as +a,.], a CR as+a,) X(s) = 
=(b+bs+...+b,s)Y(s) (27) 


definindo-se então a chamada função de transferência 
complexa (ma literatura da língua inglesa simples- 
mente função de transferência) 


(28) 


que corresponde portanto à função, no domínio 
imagem, da relação entre as transformadas de 
Laplace dos sinais de saída e de entrada do 
sistema com condições iniciais nulas. 

Deste modo conclui-se que a função de trans- 


ferência complexa é uma função racional frac- 
cionária 


bb+Hbs+...+b,s” E 


ee [4/0 SA 
à, + 2,51... a, Ss 


E (8) = 
(29-a) 


ou seja, F(s) corresponde à relação de dois 
polinómios inteiros N (s) e D (5) 


(29-b) 
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Os valores de s que anulam a função F (s) 
dizem-se zeros da função de transferência complexa 
s, (com z ==1,...,m) sendo determinados pelas 
raízes do numerador 


N (s) = bs =0 (30-a) 


ou 


N(s) =b, || (s—-s,) =0. (30-b) 


Análogamente, os polos da função de transfe- 
rência complexa são os valores de s que tornam 
a função F (s) infinita, correspondendo por 


conseguinte às raizes A (com p= 1,...,n) do 
denominador D (s) tais que 
n 
D(s) = V a, 8” =0 (31-a) 
Pp=0 
Es 
ou D(s)=a, || (6—8,)=0. (31-b) 


p= 


Em resumo, o regime de transferência do sis- 
tema é univocamente definido pelo numerador 
e denominador da função de transferência com- 
plexa, isto é 


m 
F(s) =K. ci , com K ==". (32) 
EM a 
[| (s—s,) ú 
p ==] 


Considerando o plano imagem s definido por 
um sistema de dois eixos ortogonais, sendo 
o eixo das abscissas relativo aos valores reais 
de s (ou seja, 7) e o eixo das ordenadas em refe- 
rência aos valores imaginários de s (isto é, jm), 
a representação gráfica de F(s) neste plano ima- 
gem é traduzida pelos seus polos e zeros (fig. 6-b). 

Deste modo no domínio imagem o regime de 
transferência de um sistema sem enegia é caracte- 
rizado pela representação gráfica da função de 
transferência complexa F (s) no plano imagem s. 

Mas interessa conhecer a relação entre a des- 
crição de um sistema nos domínios imagem 
e do tempo. Para se compreender essa correspon- 
dência é fundamental observar que quando o si- 
nal de entrada de um sistema sem energia for 
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a função de teste de Dirac a saida no domínio 
imagem é igual à função de transferência com- 
plexa, pois se y(t)==3(t) de (25) vê-se que 
Y(s)=1 e portanto por (27) será X(s) =F(s). 

Nestas condições a função peso representa no 
tempo o regime de transferência de um sistema 


ju Im 
Plano $ Plano Flu) 

J 0 

Cc) 


Fig. 6— Descrição gráfica de um sistema PIDy, (propor- 
cional, integral e diferencial com atraso de 1.º ordem). 
a) No tempo: por meio da função transição; 
b) No planos: através dos zeros e polos de F (s); 
c) No plano F(jw): pelo lugar geométrico 

de F (jo), ou diagrama polar. 


sem energia, pois quando se aplica à entrada 
do sistema a função de Dirac 9(t) a saída 
é a função peso g (t), e como X (s)=F(s) a defi- 
nição de transformação inversa de Laplace (24) 
conduz a LT! | X (s)) = 47! (E (s) + onde 
es (X (s)) = g (t) e portanto 

g()=<""LF(s)1. (33) 

Por outro lado a convolução é uma operação em 

que, num sistema sem energia, permite obter 
o sinal de saída x (t) relativo a um dado sinal de 
entrada y(t) através da função peso g(t), ou 
seja 

x (t) = g (t) xy (t) (34-a) 
deduzindo-se pela definição de transformação 
inversa de Laplace 


! 


niif=s | y(t-Ug(tdz. 


a 


(34-b) 


Assim a convolução traduz no domínio do 
tempo a multiplicação no domínio imagem, isto 
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é, a x(s)=-F(s).Y(s) corresponde no tempo 
x (t) = g(t) * y(t). Este resultado é altamente 
significativo na interpretação do conceito de 
função de transferência. 


3.3, Descrição no domínio complexo Z 


Aplicando-se a um sistema sem energia um 
sinal sinusoidal 


y ()=k. exp (jnt) (35) 


o sinal de saída obtém-se através da convolução 
x(t)=g(t) * k. exp (jut) (36-a) 
ou 
t 


x (t) =| g(z). kexp(jv(t— z)) dz (36-b) 


donde 


= o 


g (Z).exp(— jutz) dz— 


o 


x (t) = k. exp Goo] | 


— | g (7). exp (— jmz) a*| : (36-c) 


ÁS 


Mas para sistemas estáveis em regime esta- 
cionário as oscilações são nulas entre t e co pelo 
que o segundo intregal da expressão anterior se 
anula. Por outro lado o primeiro integral, de 
acordo com a definição de transformação de 
Laplace (23), é 


Pad 6] 


g(t). exp (—-jwnz) dz = F(jo) 


“o 


(37) 


pelo que em sistemas sem energia e estáveis 
excitados por um sinal sinusoidal y (t) se verifica 
x () = E(jo). y (8 (36-d) 


podendo-se portanto definir a resposta de frequência 
F (j») de um sistema estável sem energia 


x (t) 


+ (j 1) == 
j 19) y(b) 


(38) 


sendo F (j») a função que multiplicada por uma 
excitação sinusoidal aplicada à entrada do sistema 
reproduz exactamente o sinal de saida em regime 
estacionário. 
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A resposta de frequência F(j») equivale a 
fazer 39 == 0 na função de transferência complexa 
F (s) pois, por definição de transformação inversa 
de Laplace (24) e por (37), é £”! (g ()) = F (jm) 
e por outro lado, devido a (33), * E (8 (t)| = 
= F(s) o que se verifica apenas quando 7 =, 
ou seja s = jm, Fazer s = jm» (ou 7 ==0) corres- 
ponde a passar da transformação de Laplace 
para a transformação de Fourier, o que é possível 
se as partes reais de todos os polos da função 
F (s) forem negativas, isto é, como se estabeleceu, 
o sistema terá de ser estável. 

A resposta de frequência é uma função com- 
plexa, constituída por parte real Re e parte ima- 
ginária Im 

F (jw) = Re + j Im (39) 
quer dizer, a representação gráfica desta função 
efectua-se no plano complexo de Cauchy, obten- 
do-se o lugar geométrico de F (jw) no plano com- 
plexo (fig. 6-c), que se designa diagrama polar. 

No domínio complexo o regime de transferência 
de um sistema estável sem energia é pois carac- 
terizado pela representação da resposta de 
frequência F (j v)no plano do domínio complexo. 

Fisicamente a função F(j») equivale a uma 
variação de amplitude do sinal de entrada 


A (0) =| (jr) | (40) 
e a uma deslocação de fase 
Im CF (70) 
7 (m) == arc tg - (41) 
R (F (j 3) 


cujas representações gráficas, em escalas conve- 
nientes, reproduzem o chamado diagrama de Bode 
do sistema. 


4. DESCRIÇÃO DE SISTEMAS ATRAVÉS DO 
SEU ESTADO DINÂMICO 


4.1. Descrição no domínio do tempo 


Designa-se variável de estado u (t) de um sistema 
toda a grandeza variável interior do sistema que 
especifica o seu estado de excitação, de tal modo 
escolhida que depende apenas dos sinais de en- 
trada e ou de outras variáveis de estado. É claro 
que o processamento dinâmico de um sistema 
é independente da escolha das variáveis de estado. 
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O processamento dinâmico de um sistema 
e dos seus terminais pode ser descrito através 
das variáveis de estado por intermédio das cha- 
madas equações dinâmicas do sistema [2], que para 
múltiplas entradas e saídas se exprimem por 


ú(t)=Aut(t) + By(t 
x()=Cu(t) + Dy(t) 


(42) 


sendo A, B, Ce D matrizes coeficientes e u (t), 
y(t) ex (t) vectores relativos respectivamente às 
variáveis de estado, sinais de entrada e sinais 
de saída. 

No caso particular de sistemas apenas com 
uma entrada e uma saída as matrizes y (t), x(t), 
Be Dsão escalares y (t), x(t), bed. 

As equações dinâmicas podem representar-se 
por um esquema de blocos (fig. 7-a) onde se 
observa perfeitamente a natureza de cada uma 
das matrizes e grandezas variáveis. 


Cc) 


Fig. 7 — Esquema de blocos de um sistema, 
a) Das equações dinâmicas, 


b) Da equação de estado dum sistema linear forçado; 
c) Da equação de estado dum sistema livre. 


TECNICA N; 404 


De facto as matrizes-coeficientes das equações 
dinâmicas têm um significado bem específico: 


a) A matriz estrutural ou interna À (quadrática de 
ordem n, sendo n o número de armazenadores 
de energia do sistema) define por si só a diná- 
mica interna do sistema, isto é, o seu movimento 
próprio. Equivale ao denominador da função de 
transferência complexa. 

bh) A matriz de comando ou de entrada B (de ordem 
n.p, onde p é o número de sinais de entrada 
do sistema) indica como os sinais de entrada 
y (t) do sistema excitam. 

c) A matriz de observação ou de saída CG (de ordem 
q.n, para o número de terminais de saída q) 
caracteriza essencialmente a transferência do 
estado interno do sistema para os terminais de 
saída (pode interromper-se internamente o sis- 
tema, deixando de se observar sinais de saída, 
apenas com o anulamento da matriz O). 

d) A matriz de travessia D só não é nula se a 
correspondente função de transferência complexa 
F (s) do sistema tiver tantos zeros como polos 
(se lim F(s) = 0). Num sistema técnico (real) 

S -»+ OO 
é D=-0 pois não existe nenhum sistema que dê 
um sinal de resposta para w — “9, 
O vector das variáveis de estado 
1T 
u(lO=|u(D,u(D,..u,(D|. (43) 
pode ser interpretado como definindo uma trajec- 
tória no espaço de estado n-dimensional Rº 
(fig. 8). Neste sentido as equações dinâmicas do 
sistema representam equações vectoriais que 
indicam como o vector velocidade 


U (O) =[ (0, ú, (Dc á (| (44) 
depende do vector lugar geométrico u(t). O 
espaço de movimentos engloba ainda o eixo 
dos tempos. 

Conforme os elementos das matrizes-coeficien- 
tes assim se definem os diferentes tipos de 
sistemas: se são funções do tempo o sistema 
diz-se de tempo-não-invariante, se são constantes 
(reais ou não) trata-se de um sistema de tempo 
invariante e se são reais o sistema é linear. 
Vê-se assim que os programas de cálculo analógico 
ou digital dos sistemas lineares são essencial- 
mente simplificados. 
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No caso de um sistema não linear será 


Ú(t) = A(ut) ut) + B(ut)y(t) as) 
x (t)=C 


u() + D y (t). 


b) 


Hig.8 — Representação gráfica das equações dinâmicas 
de um sistema como uma trajectória. 
a) Sistema de 2.º ordem; 
b) Sistema de 3.º ordem. 


Num sistema livre, sem perturbação externa, 
y (t)— O para t > t,, pelo que as suas equações 
dinâmicas degeneram em 


Ú(t) = Aut) 


(46) 
xty=L wtt. 


A equação que descreve o regime dinâmico de 
um sistema, caracterizando o seu estado em cada 
instante, chama-se equação de estado do sistema 
(fig. 7-b), que no caso linear e de tempo dependente 
(não invariante) é 


ÚU()= Aut) + By(t) (47-a) 
ul) =f(u(t),y(t),t) (47-b) 


e no caso de um sistema livre (fig. 7-c) assume 
a forma simples 


Ú(t)=Au(t) (48-a) 


ou Ú(t)=f(u(t),t) (48-b) 


simplificando-se ainda mais se o sistema for 
autónomo, pois além de não ter excitações de en- 
trada é de tempo invariante, 


U (t) = f(u (t)) (49) 


sendo o sistema estacionário ou de tempo invariante 
em geral 


Úi)=f(u(t),y.t) (50) 


onde f é uma função qualquer, podendo não 
ser linear (definição de sistema não linear). 

A integração da equação de estado linear (47-a) 
pode fazer-se através da matriz de transição ou fun- 
damental 
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DP (t) = exp (A t) 51) 


e com as condições iniciais U, = u (t,), pois 
de (47-a) 
du (t) 


-— Au(t) =By(t) (47-c) 
dt 


ho) 


multiplicando à esquerda ambos os membros 
por P(t) e porque 


d du (t) 
S — At). u (t) = cam UÉ) q = 
FP (exp ( (t)) = exp ( ) E 
— hexp(— At).u (t) (52) 
e devido à comutatividade 
Â pp — A (53) 


resulta, após integração de ambos os membros 
entre t, e t e supondo as condições iniciais 
Uç ss (t,) 1 


exp(— At. u(t)-exp(—At,).u(t,) = 


a 
ii | exp(-Az).Byiz)dz (54-a) 
e multiplicando à esquerda ambos os membros 
por P (t) 


u()=u,. exp (A(t—t,)) + 


t 
+ ] exp (A(t—t,)).By(z) dz (54-b) 

adido é, 
onde a primeira parcela traduz o movimento pró- 
prio do sistema det==— — até t==t, ea segunda 
parcela corresponde à convolução das grandezas 
de entrada BY(t), representando portanto o 
movimento forçado do sistema, isto é, a varia- 
ção de estado do sistema de t, até ao instante 
de observação t devida ao sinal de entrada. 
Deste modo a transição dos estados do sistema 
é dada pela sobreposição linear dos fenómenos 
transitórios devidos às condições iniciais e à 
resposta do sistema ao sinal de excitação externa 


BY(t). 
4.2. Descrição no dominio imagem 


A aplicação da transformação de Laplace às 
equações dinâmicas de um sistema linear e de 
tempo invariante conduz às equações dinâmicas 
imagem 

sU()-u(t)=AU(s) +BY(9) 
X (5) = CU(s)+DY(s). (55) 
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Na hipótese dos sistemas sem energia, em que 
u,=u(t,)=0, pode eliminar-se o vector U (s) 
das variáveis de estado no domínio imagem, 
admitindo que existe a matriz inversa de C, pois 
da segunda equação U (s) =C =! (X(s)— DY (5)) 
e por substituição na primeira 


(sI—A) C-!(X(s)— DY(s)=BY (s). (56-a) 


Definindo-se a matriz característica imagem 


Q(s)=sI-—A, (57) 
à qual se refere o polinómio característico 
O ()=s |I-Al|,, (58) 
existe a matriz inversa 
=— 59 


desde que Q (s) seja não singular, isto é, 
se o polinómio característico não for nulo 
(det (s I-A) + 0). 

Nestas condições, multiplicando sucessivamente 
à esquerda ambos os membros da equação (56) 
por 0”! (sjeporC” !, obtém-se 


X()=(C0""'(s)B+D)Y(s) (56-b) 


ou, definindo a matriz de transferência complexa do 
sistema sem energia relativa ao seu regime diná- 
mico e dos seus terminais 


— X(g) : 
F (s) Y 65) (57-a) 
será 
F()=C0"'(s) B+D. (57-b) 


Analogamente se deduz a equação de estado ima- 
gem de um sistema linear de tempo invariante, 
com quaisquer condições iniciais, supondo que a 
matriz característica imagem é não singular 


U()=0"'(9).u, + 07! BY(s) (58) 


que se pode transformar no domínio do tempo 
u()b=u,.º. fa” (s)] + & fo” (s)BY (s)] 


(59) 


expressão que comparada com (54-b) conduz à 
conclusão que existe uma relação univoca entre 
a matriz de transição P (t) no domínio do tempo 
e a matriz característica imagem Q (s) no domínio 
imagem 
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d (= [07 "(s)) (60) 


5. CONCLUSÕES 


Pela importância do seu significado destacam-se 
os seguintes resultados expostos no decorrer 
do estudo: 


1.º Normalmente a análise dos sistemas de con- 
trole é efectuada supondo os sistemas lineares 
e de tempo-invariante. Para a compreensão das 
condições de estudo deverá ter-se presente que 
um sistema pode ser função do tempo e nestas 
circunstâncias mesmo em sistemas lineares não 
se verifica a invariância do tempo, continuando 
porém a constatar-se as leis de sobreposição e 
de homogeneidade. Nos sistemas não lineares 
por seu lado apenas se pode verificar a invariância 
do tempo (Quadro 1). 


Sistema | Linear Não linear 
Independente | Sobreposição Invariância 
do tempo Homogeneidade | do tempo 
Invariância do 

tempo 
Dependente | Sobreposição ess 
do tempo Homogeneidade 


Quadro I — Características de sistemas linear e não 
linear em função do tempo. 


2.º Existem dois processos distintos para 
descrever um sistema: através do regime de 
transferência ou por meio do seu estado diná- 
mico. Por outro lado é importante observar que 
o mesmo sistema se pode descrever em três 
domínios diferentes : tempo, imagem e complexo. 
Cada modo de descrição (Quadro II) goza de 
propriedades mais ou menos vantajosas de acordo 
com o tipo de sistema. 


3.º A equação diferencial de um sistema defi- 
ne-se pela linearização do sistema em torno do 
ponto de trabalho, traduzindo os desvios em 
relação ao ponto estacionário de trabalho, e equivale à 
identificação dos regimes de transferência dos 
terminais de entrada e de saída, pelo que a defi- 
nição no fempo de um sistema através do processamento 
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Processo de | 


descrição 
pe | ———— 
Eira is Tempo Imagem 
t s=.+jm 
Função de 
Equação transferência 
Modo de diferencial complexa 
descrição E (xx o dm X(s) 
=R(y,y' ey”) Y (s) 
| Linear Linear 
Restrições Tempo — inv. 
do sistema Sem energia 


Regime de transferência 


Estado do sistema 


Complexo Tempo Imagem 
z=Re+ijl, t 5 
Resposta Equação Matriz de 
de frequência de estado transferência 
Ecs complexa 
F (j to) ) — X(j 1») 
YG9)| U=tf(u,y,t) F (s) 
Linear Linear 


Tempo — inv. 
Sem energia 


Tempo — inv. 
| Sem energia 
| Estável 


Quadro II — Descrição de sistemas. 


de transferência refere-se sômente aos terminais e não à 
estrutura do sistema. 

4.º O regime de transferência é representado 
no tempo pela função peso g (t), no domínio 
imagem pela função de transferência complexa 
e no domínio complexo pela resposta de fre- 
quência F (j 9). 

5.º De um modo genérico a equação imagem 
de um sistema contém os valores iniciais rela- 
tivos às diferentes transformadas de Laplace da 
equação diferencial, que se devem anular para 
se poder definir a função F(s), quer dizer, a 
função de transferência complexa é apenas definível em 
sistemas sem energia, O que significa que a definição 
no domínio imagem s por intermédio do proces- 
samento de transferência só é possível se as 
condições iniciais do sistema forem nulas. 

6.º Por sua vez a descrição de sistemas no 
domínio complexo exige que os sistemas sejam 
sem energia (pois só assim se pode aplicar a con- 
volução) mas exiçe ainda a estabilidade do sistema 
(para que de facto a convolução conduza à res- 
posta de frequência). Note-se que sob o ponto de 
vista formal é costume porém definir-se por ana- 
logia, nos sistemas instáveis, numa função F (j») 
muito embora não se possa medir. Importa ainda 
observar que F(jr) não é uma função imagem mas 
sim uma função complexa, 

7.º As equações dinâmicas descrevem o comporta- 
mento dos sistemas e dos seus terminais, enquanto as 
equações de estado descrevem o comportamento dos sis- 
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temas. Com esta consideração a equação de estado 
descritiva dos sistemas tem uma forma muito 
simples, podendo simplificar-se sucessivamente 
nos casos de sistemas lineares, livres, autónomos 
e ainda de tempo invariante (Quadro III). 


Sistema Equação de estado 
Linear O = T (10, 9,:t) 
Livre ú=f(u,t) 
Tempo-inv. ú = f (u, y) 
| Autónomo u = f(u) 
Quadro HI — Equações de estado de sistemas. 


8.º A descrição no domínio imagem por meio dos 
estados do sistema não exige que o sistema seja sem 
energia, vantagem de elevado teor em relação 
à descrição através do regime de transferência 
visto que nalguns casos especiais da técnica de 
controle não é possível anular as condições 
iniciais. Analogamente a matriz de transferência com- 
plexa não exige estabilidade do sistema. 
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4 HERÓIS DA LUBR 


Audacius 


Toda a gente sabe (e se Você é Helena ainda tem 
mais razão para saber...) que foram os gregos os vence- 
dores da guerra de Tróia. 

Mas o que quase todos ignoram é que foi um chefe 
de manutenção quem, praticamente, deu a vitória aos gregos. 

Ulisses, cuja vida acabaria por ser uma verdadeira 
odisseia, nós diríamos, até, a verdadeira Odisseia, coman- 
dava o exército grego. Um exército, nesse tempo, era 
principalmente formado por infantaria, quer dizer, por tro- 
pas a pé. Daí o conhecido aforismo «debaixo dos pés se 
levantam os trabalhos». Mas Ulisses, que já andava farto 
de trabalhos (recordam-se da expressão «os trabalhos de 
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Audacius 


Hércules»? Pois vamos aplicá-la, com muito mais proprie- 
dade, a Ulisses...), teve uma ideia que só um grande ma- 
nhoso poderia segregar e acariciar: levar os troianos à 
certa metendo-lhes em casa, sem dar cavaco, um cavalo 
de pau. A ideia de Ulisses estava a caminho e acabaria por 
entrar, com um pé na História, outro na Lenda, no espólio 
cultural de todos nós. 

Mas da ideia à prática, o Cavalo nem um «salto» deu. 
Era tão pesado que não conseguia mover-se sobre as suas 
rodas de ferro. Experimentaram tirar-lhe o recheio de bem 
nutridos guerreiros gregos. Vazio, o Cavalo de Pau também 
não se movia. 

Então Ulisses consultou Audácius, o seu chefe de manu- 
tenção da época. Audácius foi categórico: era um pro- 
blema de lubrificação. «Mas o problema não se resume à 
compra de lubrificantes», avisou, de dedo engordurado no 
ar, o industrioso Audácius. «Temos de descobrir qual é o 
lubrificante indicado para cavalos de pau». 

Audácius folheou a sua papelada, que nesse tempo 
era de pergaminho, e finalmente escolheu MOBIL. 

A MOBIL não era o único fornecedor de lubrificantes 
que havia na Grécia, mas Audácius sabia que os homens 
da MOBIL eram peritos qualificados em todos os tipos de 
engenhos. 

Chegou um Delegado da MOBIL e logo iniciou os 
seus trabalhos para livrar de trabalhos Ulisses, o tal cuja 
vida foi uma verdadeira trabalheira. Começou por extrair 
de sob a sua toga um volumoso maço de pergaminhos que 
tinha por título «Manual dos Lubrificantes de Máquinas». 


Com o dedo espetado procurou no índice — «Cavalo 
de Pau». À recomendação para eixos era a Massa Lubri- 
ficante XVI — «Sovarex-Heavy Duty». 


Audácius aplicou a massa nos eixos. 

E nessa noite, com os troianos a empurrarem para 
dentro da sua própria casa o cavalo que pela calada havia 
de os escoucear, os gregos entraram em Tróia e mais uma 
página da História foi escrita. 

Hoje em dia, os agentes de compras como Audácius, 
procuram não só bons produtos, mas também boa assistência. 

Por exemplo: o «Manual dos Fabricantes de Máquinas» 
é o guia de que os Delegados da Mobil se servem para 
recomendarem o lubrificante mais adequado a cada órgão 
de uma máquina. Este valioso livro contém uma relação 
dos mais importantes fabricantes de máquinas industriais 
e descreve os equipamentos por eles fabricados. Para cada 
parte móvel dessas máquinas há uma recomendação espe- 
cial de lubrificacão. 

Não é preciso um sujeito chamar-se Audácius para 
saber como orientar-se na recomendação dos lubrificantes 
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TRATAMENTO DE RESÍDUOS 


RESUMO 


O aumento da produção de resíduos obriga a estu lar 
novos métodos de tratamento e a melhorar os existentes, 

Este artigo descreve o estado actual do problema e estuda 
alguns processos sub o ponto de vista econômico, 


1 — INTRODUÇÃO 


O tratamento da água poluída dá origem a 
resíduos que é necessário armazenar em pura 
perda de espaço, ou transportar para locais mais 
espaçosos, como o oceano. 

O crescimento das quantidades de resíduos 
formadas num dado intervalo de tempo aumenta 
extraordinariamente à medida que aumentam 
o nível de vida e a industrialização, necessitan- 
do-se de novos métodos para que a área dispo- 
nível não seja em pouco tempo exígua. 

O transporte destes resíduos é caro e por isso 
a hipótese de redução de volume tem-se apresen- 
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SYNOPSIS 


Rapid depletion of disposal areas makes imperotive to 
dvelop new approaches for handling wastes and to improve 
established methods, 

This paper looks for the state-of-the-wrt and studies some 
motels under the economic point of vie, 


tado como a mais tentadora para a resolução do 
problema do armazenamento de resíduos. 

No entanto, à maior redução de volume, 
corresponde naturalmente um maior custo, pelo 
que este valor terá que ser compensado pelos 
custos de transporte e armazenamento. 

Como exemplo observe-se o diagrama da 
figura 1 preparado pela firma americana na 
Syska & Hennessy para um hospital com 
600 camas produzindo por dia cerca de 14 tone- 
ladas de resíduos sólidos correspondentes a 
190 m” | dia. 

Nela estão indicados os vários processos consi- 
derados. A escolha do processo mais económico 


EVOLUÇÃO DOS CUSTOS INICIAIS E OPERATÓRIOS COM O AUMENTO DA REDUÇÃO DE VOLUME 
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